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"^ ■ Abstract 

The Kramers turnover problem, that is obtaining a uniform expression for the rate of escape 

Oj, of a particle over a barrier for any value of the external friction was solved in the eighties. Two 

formulations were given, one by Mel'nikov and Meshkov (MM) (J. Chem. Phys. 85, 1018 (1986)), 

,i^ ■ which was based on a perturbation expansion for the motion of the particle in the presence of 

O 
cy^ ! friction. The other, by Pohak, Grabert and Hanggi (PGH) (J. Chem. Phys. 91, 4073 (1989)), 

O . 

'c/3 ' valid also for memory friction, was based on a perturbation expansion for the motion along the 

>.' 
r^ ' collective unstable normal mode of the particle. Both theories did not take into account the 



temperature dependence of the average energy loss to the bath. Increasing the bath temperature 



^ I will reduce the average energy loss. In this paper, we analyse this effect, using a novel perturbation 

Cf^ ' theory. We find that within the MM approach, the thermal energy gained from the bath diverges, 

the average energy gain becomes infinite, implying an essential failure of the theory. Within the 
f~^ i PGH approach increasing the bath temperature reduces the average energy loss but only by a 

m ■ 

finite small amount, of the order of the inverse of the reduced barrier height. This then does not 
seriously affect the theory. Analysis and application for a cubic potential and Ohmic friction are 
^ , presented. 
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I. INTRODUCTION 

One of the outstanding problems of the second half of the twentieth century in rate 
theory was the solution for the rate of escape of a particle trapped in a potential well over 
an adjacent barrier whose motion is subject to a Gaussian random force and associated 
frictional force. In 1940, Kramers [l| solved this problem in two limits. In the underdamped 
limit, when the friction is weak he found that the rate increases linearly with increasing 
friction. In the opposite strong damping limit, the rate decreases inversely with the friction 
strength. Kramers however did not manage to derive an expression for the rate which is valid 
for the whole range of friction. This problem is known as the Kramers turnover problem. 

An important milestone in solution of this problem was given by Mel'nikov and Meshkov 
(MM) in 1986 J2|, |3|. They introduced into the problem the concept of a conditional proba- 
bility for the particle initiated at the barrier with energy E to return to it with energy E'. 
Using perturbation theory with the friction strength as the small parameter they derived an 
explicit Gaussian expression for the kernel, in which the central quantity was the average 
energy lost by the particle to the bath as it traversed from the barrier to the well and back. 
This energy loss was temperature independent. Given the kernel, they wrote down a master 
equation for the energy flow and upon solving it were able to derive an expression for the 
depopulation factor rate which went smoothly from the weak damping limit where it was 
linear in the friction strength to unity for strong damping. They then wrote the full rate 
expression as a product of three terms, the "standard" transition state theory expression 
for the rate, the depopulation factor, and Kramers' spatial diffusion factor which went from 
unity at weak damping to being inversely proportional to the damping when the friction 
was strong. 

Kramers' original theory dealt with Ohmic friction, that is the white noise case. Grote 
and Hynes noted in the early eighties that in Chemistry one should expect the environment 
to have memory, since the typical time scales of a molecular surrounding are similar to that 
of the reacting system. They then proceeded to solve the Kramers rate problem in the 
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presence of memory friction, in the spatial diffusion limited regime j^, that is in the regime 
of moderate to strong friction. Carmeli and Nitzan [Sj proceeded to do the same for the 
underdamped regime. However the turnover problem in the presence of memory friction 
remained unsolved, even after the seminal result of Mel'nikov and Meshkov. 



The solution for the turnover problem, without the need for the ansatz that the rate is 
a product of three terms, was derived in a few steps. First, PoUak 6| realised in 1986 that 
Kramers' problem in the spatial diffusion limit may be derived using variational transition 
state theory, based on a Hamiltonian description of the dynamics, in which the system is 
bilinearly coupled to a harmonic bath. In the vicinity of the barrier, the Hamiltonian is a 
quadratic form and so can be diagonalized. The normal modes then include an unstable 
collective mode and stable bath modes. The unstable mode is a linear combination of the 
system coordinate and the original (coupled) bath modes. Then Grabert [TJ, using the 
discretized oscillator model derived a generalization of Mel'nikov and Meshkov's turnover 
theory by considering the dynamics of the unstable normal mode, instead of the system 
coordinate. The continuum limit version of this theory was then derived in Ref. 8| and is 
known as PGH theory. Compared to Mel'nikov and Meshkov's theory, PGH theory has two 
advantages: a) it is derived, there is no ansatz in the derivation; b) The result is also valid 
for memory friction. Practically though, the expression for the energy loss of the particle 
appearing in Mel'nikov and Meshkov's theory is much simpler to implement. 

Both theories ignored the fact that the average energy loss should depend on the bath 
temperature. As the bath temperature increases, the fluctuations increase and the bath not 
only dissipates energy from the system due to the frictional force, it also pumps energy back 
into the system. As long as the energy of the system is much higher than the thermal energy 
{ksT) one expects that the particle will on the average lose energy to the bath, however as 
the temperature is increased the average energy lost will be reduced. This is the topic of this 
present paper. We will use a perturbation expansion to derive a term for the energy gained 
by the particle during one traversal, using the PGH approach, based on the motion of the 
unstable normal mode and the MM approach, based on the motion of the system coordinate. 
We will show that for the MM theory the average energy gain diverges, since the particle 
feels the fluctuations of the bath during all of its traversal. This is then an essential failure 
of the MM theory. For PGH theory though, to leading order in the perturbation theory 
the average energy loss is a sum of two parts. One is the temperature independent energy 
loss, due to the frictional force. The other is a fluctuation induced energy gain which grows 
linearly with the temperature of the bath. This energy gain is finite, since in the vicinity 
of the barrier the unstable mode motion decouples from the bath. Inclusion of the energy 
gain term in the depopulation factor is then shown to lead to a correction which is inversely 



proportional to the reduced barrier height, and so typically, can be neglected. 

In Section II we consider the average energy loss and energy gain within PGH theory. 
Then in Section III we consider the average energy loss and gain within MM theory and 
show that for Ohmic friction, the average energy gain, at any finite temperature, diverges. 
In Section IV we analyse the effect of the energy gain on the depopulation factor within PGH 
theory. A specific analytic example for a cubic oscillator and Ohmic friction is considered 
in Section V. We end with a discussion of implications and further improvements of PGH 
theory. 

II. THE ENERGY LOSS IN PGH THEORY 

A. Preliminaries 

The classical dynamics of the generic system is that of a particle with mass M and 

coordinate q whose classical equation of motion is a Generalized Langevin Equation (GLE) 

of the form: 

Mq + ^^^ + M f dt''y {t - t') q {t') = F (t) . (2.1) 

dq Jo 

F (t) is a Gaussian random force with zero mean and correlation function 

{F{t)F{t')) = MkBT-f{t-t'). (2.2) 

7 (t) is the friction function, fc^ is Boltzmann's constant and T is the temperature. The 
potential is assumed to have a well at g^ with frequency Ua and a barrier at g = which 
separates the well from either another well or a continuum. The harmonic (imaginary) 
frequency at the barrier top is denoted as u^. The potential is then written as 

V{q) = -^Mu^'q' + V,{q) (2.3) 

and Vi (q) is termed the nonlinear part of the potential function. Kramers' turnover theory 
provides an expression for the rate of escape of the particle from the well to the adjacent 
well or to the adjacent continuum as a function of the friction function, the temperature, the 
barrier height V^ = V [0) — V [qa), the mass of the particle and properties of the potential. 



Is 



It is well known J9| that the GLE may be derived as the continuum limit of the equations 
of motion of the particle bilinearly coupled to a bath of harmonic oscillators, with mass 



weighted coordinates Xj and momenta pj-. 



2 1 ^ 



j=i 



p^ + W^ I Xi 



Mc 



u 



.^1 



(2.4) 



where p is the momentum of the particle. The friction function is identified as 

7 (^) = y] ^ cos (ujt) 

j = l J 

and the random force is a function of the initial conditions of the bath oscillators: 
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Here, the subscript denotes the initial time. One readily finds that averaging the random 
force with the canonical distribution exp(— /3/7) with (/3 = 1/ (fc^T)) implies that it is 
Gaussian with zero mean and correlation function as given in Eq. 12. 2[ 

Around the barrier top, the Hamiltonian has a quadratic form and so may be diagonalized. 
The details of this normal mode transformation may be found for example in Ref. 10|]. The 
Hamiltonian is written in terms of the system (unstable) mass weighted normal mode p and 
momentum Pp and stable bath normal mode coordinates i/j and momenta Py, as: 



H 
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N 



U'p' + VM + \Y.Vl+X]y^ 



i=i 



(2.7) 



with Aj, the frequency of the j-th normal mode. The system coordinate q is expressed in 
terms of the normal modes as 



N 



Mq = uoop + ^ Ujoyj 



(2.8) 



where Ujo is the projection of the system coordinate on the j-th normal mode. The nonlinear 
part of the potential Vi (g) couples the motion of the unstable normal mode (p) motion to 
that of the stable normal modes {Vj). The unstable normal mode barrier frequency A-^ 
is expressed in terms of the system barrier frequency and Laplace transform of the time 
dependent friction {■y (s)) via the Kramers-Grote-Hynes relation [l|, |j]: 



A^^ + 7 (A^) A^ 



u 



t2 



(2.9) 



The matrix element mqo for the projection of the system coordinate on the unstable mode is 
also expressed in terms of the Laplace transform of the time dependent friction |8|: 



■^00 



1 l{X^) ^djis) 



(2.10) 



These last two relations imply that the unstable mode barrier frequency and the projection 
amplitude mqo are known in the continuum limit. 

The assumption of weak coupling between the system and the bath is expressed as |8|: 

uI = 1-uIq<^1 (2.11) 

which equivalently implies that the normal mode transformation elements Ujo,j ^ are 
small. 

The equation of motion for each of the stable normal modes is: 

^. + % = -^K (^^ (^nooP + E-.oy,)) (2.12) 

while the equation of motion for the unstable normal mode is: 

p - A^V = -^V; (-^ Lop + f: u,oy] j . (2.13) 

B. Perturbation theory 



The motion for each of the stable bath modes may be represented in terms of a power 
series in the small parameter Ujo- For our purposes it suffices to consider the first two 
terms in such a series, so that for example the coordinate of the j-th stable mode has two 
components |/jo,t and yji^t of order m°q and mJq respectively. Similarly the system unstable 
mode motion may be expanded in a perturbation series such that the first two terms for the 
coordinate of the unstable mode po,t and pi^t are of the order m° and u\ respectively. The 
equations of motion for the components of each of the stable modes is 

%o,t + A>,o,i = (2.14) 

ftu + Vft„ = -^n(^ (2.15) 



Similarly the equations of motion for the first two components of the unstable normal mode 
are 

Po,t 

and 



-\''pot = -^Vl(^^^] (2.16) 



Pi,t - A*Vi,t = -^^i" ( "77^ ) ( "ooPi,t + Xl^ioyio,* j • (2.17) 

A central quantity in Kramers' turnover theory is the energy gained by the bath as the 
unstable mode motion moves from the barrier region to the well and back. The stable mode 

bath energy is by definition 

1 ^ 






The change in time of the bath energy including all terms up to u'^q is readily seen to be: 



dE, ^ 



dt 



E 



Ujo{yjo,t + yji 









+0{u%) (2.19) 

so that to this order it suffices to consider only the zero-th and first order solutions for both 
the stable modes and the unstable mode. 

The solution for the zero-th order equation of motion for the j-th stable mode is that of 
a harmonic oscillator 

%-,o (t) = yj,o (to) cos [A, (t - to)] + ^^^ sin [A, (t - to)] . (2.20) 

It is useful to define a Gaussian random "noise function" $ (t) as: 

TV 

$(t) = ^M,o%o.*- (2-21) 

i=i 

When averaging over the uncoupled bath Hamiltonian using the thermal distribution 
exp {—^Eb (to)) one readily notes that the mean of $ (t) vanishes and the noise correla- 
tion function is 

($ (t) $ it')) = ^Y.iT^os [A, (t - t')] ^-^Kit- t') (2.22) 



and this defines the "normal mode friction kernel" Kit). Using properties of the normal 

n 

mode transformation (see for example Eq. 2.17 of Ref. [10[) one may readily express the 
Laplace transform of the kernel as 

j>(A-( ^^00 I g + 7(g) \ .2 2S^ 

showing that it is defined in terms of the Laplace transform of the friction function and thus 
in the continuum limit. 

The solution for the first order equation of the j-th stable mode is (taking tg -^ — oo): 

_u^ r ^,.Bin[A,(t-tO] fu^^\ 



y M J -oo 'V 

and it is independent of the noise $. 

The zero-th order equation of motion for the unstable mode is that of a conservative 
system with a nonlinear potential. In some cases, the trajectory for the motion from the 
barrier top to the well and back to the barrier top is known analytically. In any case, the 
period of this orbit is infinite, it starts in the infinite past at the barrier and returns to 
it in the infinite future. Henceforth we will limit all solutions to these conditions, so that 
for example the initial time is taken as — oo. The solution for the first order equation of 
motion for the unstable mode is somewhat more complex, however readily solved by using 
the conservation of the overall energy of the system and the bath, as shown in Appendix A. 
One finds that it is linear in the noise: 

ft' 

Po,t' 

For later purposes, we note that the thermal average of the product of the first order solution 
with the time derivative of the noise is: 






,„i(()\ = -^ /' M^ f df'-^zn^v; h"p^"y] . (2.26) 

"■ ^V VMl3j-^ Plt'J-o. dt dt" '\ -/M ' 



C. The thermally averaged energy loss and its variance 

The thermally averaged energy loss is obtained by averaging the energy change over the 
thermally distributed initial conditions of the bath exp(— /3_Eb (— oo)). This implies that 



the thermal averages of the stable mode coordinates and velocities vanish. The coordinate 

correlation functions are: 

cos [A,- (t — t')] ^ 

{y,o,tyko,t') = — ^-j^2 — -^j'' (2-27) 

where 6jk is the Kronecker " delta" function and from this relation one obtains by differenti- 
ation with respect to the time the velocity correlation functions. This implies for example, 
that /<i>(t)<l>(t)\ =0. 

The average energy loss of the system, defined as the gain in the energy of the bath as 
the unstable mode motion goes through one cycle, starting at the barrier and returning to 
it, is obtained by integration of Eq. 12. 191 over the unperturbed barrier top trajectory for the 
unstable mode, that is the time goes from — oo to oo. It is then readily found to be 

(AE) = r dt (^) = DEi - DE2. (2.28) 

Using the first order solution for the stable mode motion (Eq. I2.24p one finds that the 
systematic component of the energy loss DEi is: 

2M J_^ J_^ \ s/M J Stat' ' V Vm y 

and this is the term that also appears in PGH theory. The new result is that to the same 
order in the perturbation strength, there is also a temperature dependent contribution to 
the energy loss: 

M2/3 J_^ J_^ ' \ ^/M J ^°'* dt ' \ y/M J ^°'* Jt" pI 



Po,t' 
(2.30) 

which expresses the fact that when the bath is heated it may also transfer energy to the 

system. It is this temperature dependent reheating of the system which is missing in PGH 

theory. 

The variance of the energy loss is by definition 

a^ = {{AE - {AE)f) . (2.31) 

To leading order we use only those terms which are up to quadratic in Ujq. This implies 
that the term (AE) which is fourth order in Ujo may be ignored. The leading order term 
in the energy loss is then identical to that found in PGH theory, that is: 

a^ = ^DEi. (2.32) 



III. THE ENERGY LOSS FOR THE SYSTEM COORDINATE 

The turnover theory of Mel'nikov and Meshkov was based on a perturbation expansion 
in which one considers directly the motion of the system coordinate. The small parameter 
is taken to be the noise F (t) and we note that the friction function is second order in the 
noise. The system motion is separated into three parts 



qt = qo,t + qi,t + q2,t + O (F^) 



(3.1) 



which are of the order F , F and F respectively. They obey the equations of motion: 

= Mqo,t + V'iqo,t), (3.2) 

F{t) = Mq,^t + r'iqo,t)qi,t (3.3) 

= Mq2,t + V" iqo,t) q2,t + -V^'^ iqo,t) qlt + M / dt'^ (t - t') qo,t'. (3.4) 

The time derivative of the kinetic energy of the system, including all terms up to order F^ 
is then: 



dE 



qo,t 



-V iqo,t) - V" (go,*) {qi,t + q2,t) - V^^) (go,*) q^ -M [ dt'^ (t - f) %,, + F (t) 

^ J-to 



Hu [-V {qo,t) - V" (go,i) gi,i + F (t)] - qt,2V' {qo,t) + O {F') 



(3.5) 



The energy loss is obtained by integrating using the zero-th order trajectory that is initiated 
at the barrier and returns to the barrier. Using the perturbative equations of motion given 
in Eqs. I3.2ti3.4l one finds that to order F"^ the energy loss depends only on the zero-th and 
first order contributions: 



-AF„ 



~ , dE 
dt^ 
dt 



— oo 
oo 



/oo nt 

dt -Mqo,t / dt'-i (t - t') go,i. + (go,* + gi,t) F (t) 
oo L J —oo 

Averaging over the random force then gives two contributions: 



O (F^) . (3.6) 



- (A^.) 



dt 



Mqo,t / rft'7 (t - t') qo,t' + (gi,tF (t)) 



+ O (F^) . 



(3.7) 



The first contribution on the right hand side, which is independent of temperature, is the 
same as derived in the turnover theory of Mel'nikov and Meshkov. The second contribution 
is though missing. 
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As in the PGH approach, one may derive an exphcit expression also for the second 
temperature dependent contribution to the energy loss. This is readily carried out using 
the Hamiltonian representation for the Langevin equation, as given in Eq. 12.41 For this 
purpose, the motion of each bath mode is also treated with perturbation theory, where the 
small parameter is the coupling constant Cj. The zero-th and first order contributions to the 
motion of the j-th bath mode are then denoted as Xjo,t and Xji^t respectively. The zero-th 
order motion is that of the uncoupled j-th bath oscillator. The first order correction to the 
j-th bath oscillator motion is readily found to be: 



xji,t = VMcj / dt' '-^ ^go.t'. (3.^ 



Following the same derivation as in Appendix A, one notes that energy conservation 
implies that to first order in the coupling coefficients Cj-. 

N 

= -^^J^ + ^' ('?o) ?i + X^ [PjoPji + (^]xjoXji] - Fqo. (3.9) 

i=i 

This then gives a first order in time equation of motion for the first order correction to the 

system coordinate which is readily solved: 

qi,t = qo,tf dt'-^ [ dt"qo,t"Fit") (3-10) 

We thus find that: 

»t 1 rt' 

-oo ^^^%,t' J-oo 

.. , dt'qo,t'l{t-t'). (3.11) 

P(l0,t J-oo 

The energy gained by the system from the bath is thus: 



. M din go. f ,. 1 f .„. ., .,. 

h,tF{t)) = — — -— go,t / dtj—^ / dt qo,t"l{t-t 
/3 dt J_^ Mq^,, J_^ 



1 >■' 



/oo 
dt{qt,iF{t)) 
■oo 



oo 

J\^Jf /»00 /"OO i*t -I 



dtqo,t / dt"qo,t"l{t-t") / dt' 

P J-oo J-oo Jt" ^^^Qo.t' 



0,t' 



oo -I rt 



+ / "^^IT— dt'%,t'l{t-t'). (3.12) 

J-oo P<l0,t J-oo 

Consider then the case of Ohmic friction, that is: 

7 (t - t') = 27(5 (t - t') (3.13) 
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where S (x) is the Dirac "delta" function. The first term on the right hand side of the second 
equahty in Eq. 13. 121 vanishes but the second term diverges: 



(AE,), 



7 



dt — )• oo. 



(3.14) 



'2,Ohmic a 

This is then an essential failure of the Mel'nikov and Meshkov solution for the turnover 
problem. Thermal fluctuations of the bath give an infinite amount of energy to the system. 
The same divergence does not occur in PGH theory where one considers the motion of the 
unstable normal mode. In the vicinity of the barrier, the unstable mode is decoupled from 
the bath and the bath can no longer affect the motion. In the Mel'nikov Meshkov theory, 
even arbitrarily close to the barrier top, the system coordinate remains coupled to the bath. 
Due to the infinite time duration of the motion, the bath can then provide an infinite amount 
of energy and this causes the divergence. 



IV. PGH TURNOVER THEORY 



A. The energy loss conditional probability 



A central object in the turnover theory is the probability kernel P {E'\E) that the system 
originates at the barrier with energy E and returns to it with energy E' . From Eq. 12. 191 one 
readily finds that the change in energy may be written as: 



AE = E-E' = DEi + / dtWi 



$(t) 



— oo J — oo 



+ / dt dt'W' 



$(t),<l>(t') 



dtW', 



$ (t) , $ (t) 



with 



Wi 



<j>(t) 



1 y;^ooM^)')^(,) 



W. 



w. 



-$ (t) - 

M ^' dt 



$ (t) , $ (0 = -$ (t) ^ k f^") 1 $ (0 ^v[ f^^^ 

^^ ^ ^ .. W ^, 1 \^i)\ ^ ' dt' '\ Vm 



$ (t) , $ (t) 






(4.1) 



dt"^ 
t' Po,t 



(4.2) 

(4.3) 
(4.4) 
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The thermally averaged probability P {E'\E) that the particle ends with energy E' after 
being initiated with energy E is then by definition 



P{E'\E) = 
sIe' -E + DEi + 



dt 



Wi 



<i.(t) 



w. 



$ (t) , $ (t) 



dt'Wo 



$(t),<l>(t') 



/3 

(4.5) 



where the average is over the distribution 

TV 



7 = 1 ^ ^ 



(4.6) 



with py. = ijj and we used the shortened notation yjQ = i/jq (tg), etc.. 

If one considers only the first two terms in the energy loss (ignoring the terms with W2 
and W3) one regains after performing the thermal averaging, the PGH kernel 



Po (E'lE) 



/3 



1/2 



exp 



/3{E' -E + DEi) 



(4.7) 



A-kDEiJ ' \ 4DEi 

One can also explicitly perform the averaging over all modes with the full expression (in- 
cluding the terms with 14^2 and W3), since it involves only Gaussian integrations however 
we have not found a way to express the final result in the continuum limit. Instead we will 
consider the lowest order perturbation theory limit. 

For this purpose we note that the conditional probability, to lowest order has to obey 
three conditions. The first is normalization, that is 

/oo 
dE'P{E'\E) = l. (4.8) 

-00 

The second is that the average change in energy must be given by the average energy change 

found in the previous section, that is 

/oo 
dE'P {E'\E) {E - E') = DEi - DE2. (4.9) 

-00 

Thirdly, the conditional probability must obey detailed balance, that is 



P {E'\E) exp {-/3E) = P {E\E') exp {-/3E') . 



To simplify, we use henceforth the following reduced notations: 



e = I3E, 5i = ^DE,, 82 = PDE2, y = e' - e. 



(4.10) 



(4.11) 
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The leading order correction to the conditional probability, which also obeys detailed balance 
must be of the order y^. We then write the conditional probability as: 

P{e'\e) = a{l-hy^)Po{e'\e). (4.12) 

The two coefficients are determined by the normalization condition and the known thermally 
averaged energy loss: 

1/2 



l^£.,a(l-.,^)(J^) e.p(-(l4^) (4.13) 

..-.. ^ -jyyai,-,y^)[^)"\Jjl±^). (4.14) 



One then finds that: 



P {e'\e) = [l + ^^{2 + 6,~y^y^ P,{e'\e) . (4.15) 



B. The depopulation factor 



The depopulation factor given in the turnover theory is 



-^ /-oo In 



l-P(r-i) 



where -P (''" ~ f ) is the Fourier transform of the conditional probability kernel. One finds: 
P(r--j = / dy e^^ \i{T ~ -)y\ P {y) 

l + 52(r^ + i))exp(-5i(r^ + i)). (4.17) 

Assuming that the correction is small, the depopulation factor is then expanded to 

T^Toexp(-|^(l-erf(^))) (4.18) 

where Tq is the depopulation factor using the zero-th order kernel in Eq. 14.161 One notes 
that the thermal correction 82 reduces the depopulation factor and thus the rate. The 
reduced energy loss of the particle, effectively makes it more difficult to energize the particle 
and the rate is reduced. As we shall see considering a specific example below, the ratio j- is 
typically of the order of 1/{PV'^) so that it is small and does not affect the barrier crossing 
rate appreciably. 
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V. AN EXAMPLE - OHMIC FRICTION AND ESCAPE FROM A CUBIC PO- 
TENTIAL WELL 



For Ohmic friction, the memory kernel in the normal mode representation as taken from 
Ref. lul is: 



^ (t) = ^ (exp i-XH) - ^ exp (-Alt) + 1 + ^) " ;^ (5-1) 



with 



2 \ 1/2 



For a cubic potential 



V{q) 



Mu^' 



one finds 






(5.2) 



(5.3) 



11| that the solution for the zero-th order equation of motion for the unstable 
mode (Eq. 12161) is: 



Po,t 



U^OqU 



t^oshM^t 



SO that 



BE2 = -7 dx dAx (1 - exp (-2Ax-^jj exp {-2Ax 

sinh (x + ^) sinh [x — ^) 



cosh^ (x + ^) cosh^ (^ - ^) 
where we used the notation 



M(x + — j -M ix- — 



-* , cosh^x _/A^ \ , /A^ „. 

rfx 2— = M I — t 1 - M I — t" ) . 

Ai+// sinh X 



In the limit of weak friction ^ ^ 1 this reduces to 



DE2 (7 ^ 0) 



57 



In this limit, the temperature independent part of the average energy loss is [11 1 

3671/^ 



DEi (7 -^ 0) 

so that the ratio of the two terms is 

DE2 (7 ^ 0) 



5 ut 



25 



DEi (7 -> 0) 36/31^^ 



(5.4) 



(5.5) 



(5.6) 



(5.7) 



(5.8) 



(5.9) 
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which is a small number. As noted above, the temperature induced correction may be 
ignored in this limit. 

It is also of interest to estimate the temperature dependent energy gain in the opposite 



limit that is when ^ ~ ^ ^ 1. One readily finds that 



DE2 (7 -)■ 00) 



/?■ 



In this same limit the temperature independent energy loss is llj 

27 



so that: 



DEi (7 -^ 00) 
DE2 (7 -;■ 00) 



35 

175 



(5.10) 



(5.11) 



(5.12) 



DEi (7 -> 00) 27/3V^^ 
and this is not necessarily small. Strictly speaking, in this limit the perturbation expansion 

is invalid, however the result does point out that when the friction is not too weak one may 

find discrepancies between the "standard" PGH theory energy loss and the true energy loss, 

which originate from the fact that the thermal bath can also contribute to decreasing the 

energy loss of the system. 

0.25 



0.2 - 



< 0.15 |- 



< 



0.1 - 



0.05 - 



T 



T 



T 



_L 



0.0 



0.05 



0.1 



0.15 



0.2 



FIG. 1: Frictional dependence of the ratio of the thermally induced energy gain ()5.5p and the energy 
loss (j2.29p for increasing barrier height (top to bottom): PV^= 4 (dotted), 6 (short-dashed), 8 
(dashed), 16 (solid). 

The above results are illustrated in Figs. [H |2] where we show the thermally induced 
energy gain relative to the energy loss and the reduction of the rate due to a reduction 
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of the depopulation factor (14.161) compared to the PGH result. According to the above 
discussion, in the domain of low friction the ratio DE2/DE1 is mainly determined by (3V^ 
with the tendency to be more sensitive to dissipation for lower (reduced) barrier heights. 
The impact of the energy gain on the depopulation factor decreases with growing friction 
and increasing barrier height as indicated in (I4.18p : for fixed barrier height (52/(5i is almost 
constant while 61 ~ ■yV^ grows so that T/Tq — ;■ 1. In typical experimental situations with 
/3V^ of the order of 10, small deviations from the PGH prediction may be found in the weak 
damping regime 7/w < 0.1. 
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FIG. 2: Ratio of the depopulation factors T/Tq according to (|4.16p with (I4.15P and without (14.71 
the thermal correction to the energy transfer kernel. 



VI. DISCUSSION 



The solution of the Kramers turnover problem in the 1980s has been a cornerstone result 



of thermal rate theory 12|, ll3| . The rate expression obtained within PGH theory has been of 
practical use in a variety of fields such as chemical reactions [M, [l5| and Josephson junction 
physics [l6|, |l7|. However, it is incomplete as we have shown in this work since it does 
not include the fact that the temperature of the surrounding bath degrees of freedom will 
typically reduce the average energy lost by the system, due to the thermal fluctuations. 
Within a consistent perturbative treatment, this "heating" contribution should be taken 
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into account on the same footing as the dissipative part. One then finds that the thermal 
fluctuations indeed lead to a reduction of the energy loss that a particle experiences when 
traversing the well of a barrier potential. In contrast, the variance is affected only in higher 
order of perturbation theory. Consequently, the depopulation factor of the rate is reduced 
compared to the PGH prediction. 

This result cannot be obtained within the Mel'nikov and Meshkov formulation where the 
heating contribution diverges. While the corresponding depopulation factor has often been 
used in practical calculations due to its simpler structure, our analysis demonstrates the 
inconsistency of the approach. Discrepancies between the corrected rate expression of the 
PGH theory and the original one are expected to be detectable in the low friction regime 
and in systems with moderate (reduced) barrier heights. 

In this paper we have limited ourselves to the lowest order consistent contribution to 
the rate expression. Within the PGH framework, it is possible to go to higher order and 
thus obtain a further temperature dependent consistent correction also to the variance of 
the energy transfer. One may use the same perturbation theory to also obtain improved 



dynamical estimates to the rate, using the formalism of Ref. 18|. 
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Appendix: Unstable mode motion in first order 

In this appendix we note that by employing energy conservation of the composite system 
and bath one obtains a flrst order in time equation of motion for the flrst order correction 
to the unstable mode motion. Using the representation of the Hamiltonian in terms of the 
normal modes (Eq. 12. 7p one readily flnds that to flrst order 

Ppi,tPpo,t = A^ po,tpi,t j= V^ (^W) ~ ^ [Pyj^,tPy,o,t + ^jVjhtyjoA ■ (A.l) 

V \ V / j—i 
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Noting that 

Pi,t = Ppi,t (A. 2) 

we obtain a first order in time equation of motion 

Pht = Yf P'^^o.t] Pi,t - T=Vi -y=- - — 2^ [Pyn,tPy,o.t + \l/ii,t%o,tJ • (A.3) 

at po,tvM V V M / Po,t ~^ 



The solution for the first order correction to the motion of the stable modes has been given 
in Eq. 12.241 Noting that at the initial time to the potential term Vi i "°'^'* j vanishes, one 



finds that the solution for the first order term for the unstable mode motion is: 

Pl,t 



^^[\,J^f\f^jn^y^hj^\ (A,4) 

i,. I'le Ju, VM dt" ^\ VM J 



and it is linear in the noise $. It is an instructive exercise to show by direct differentiation 
with respect to the time that this solution is also a solution of the second order equation of 
motion, as given in Eq. 12.171 

One may follow this procedure in principle to also obtain the higher order solutions. 
Specifically, given pi ^ one obtains an explicit result for the second order contribution to the 
j-th stable mode yj2,t- Then one uses again the conservation of energy to obtain a first order 
in time equation for p2,t which is readily solved, and so on. 
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